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Using the canonical Monte Carlo simulation technique, we study a Regge calculus model 
on triangulated spherical surfaces. The discrete model is statistical mechanically defined 
with the variables X, g and p, which denote the surface position in R, 3 , the metric on a 
two-dimensional surface M and the surface density of M, respectively. The metric g is 
defined only by using the deficit angle of the triangles in M. This is in sharp contrast 
to the conventional Regge calculus model, where g depends only on the edge length of 
the triangles. We find that the discrete model in this paper undergoes a phase transition 
between the smooth spherical phase at b— too and the crumpled phase at b— >0, where b is 
the bending rigidity. The transition is of first-order and identified with the one observed 
in the conventional model without the variables g and p. This implies that the shape 
transformation transition is not influenced by the metric degrees of freedom. It is also 
found that the model undergoes a continuous transition of in-plane deformation. This 
continuous transition is reflected in almost discontinuous changes of the surface area of 
M and that of X(M), where the surface area of M is conjugate to the density variable 
P- 
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1. Introduction 

A surface model for membranes is defined by a mapping X from a two-dimensional 
surface M to R 3 l l l ^PK I Qn iy ^ e ma ppi n g x but also the metric g of M is 
assumed as dynamical variables of the model in its statistical mechanical study ^1. 
However, the variable g is always assumed to be the Euclidean metric g a b = Sab or 
the induced metric g a b = d a X^dbX^ of the mapping X in the numerical studies 
that have been performed so far In those studies, the variable g is fixed in the 
partition function, and as a consequence a role of g in the phase structure of the 
model remains unclear. 

The metric g of M can be discretized as a y ariable by using the edge length of 
triangles in the Regge calculus approach to the triangulated surface model. 
The model is defined by the Hamiltonian S — Si + frS^, where Si and S2 are the 
Gaussian bond potential and the bending energy, and b[kT] is the bending rigidity 
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| iu | ii | iz | j n a pp r oach, one can see whether or not g influences the discontinuous 
transition between the smooth phase at b — > oo and the collapsed phase at 6^0. 
This transition and related phe nomena have long been studied numerically as well 
as analytically by many g r o U p S l l3 | 14 | 15 | 16 | 17 | 18 | 19 | 2() | 21 | 22 | 23 | 24 | 25 | 26 | 27 | 

Such a Regge calculus model was recently studied in Ref. [551 where M is em- 
bedded in R 2 and the metric g a b is slightly extended from the conventional Regge 
metric by incorporating a degree of freedom for a deficit angle of the triangles in 
M. One of the non-trivial results reported in Ref. [28] is that the deficit angle varies 
almost discontinuously at the transition point, although the discontinuity is very 
small compared to the deficit angle itself. 

However, it is still unclear whether or not such a discontinuity in the internal 
geometric quantities can also be seen in the case where the edge length of triangles 
is fixed while the deficit angle is varied. Moreover, it is also nontrivial whether 
the fluctuating metric influences the phase structure of the model with a constant 
metric. 

In this paper, we study a surface model with a metric variable, which depends 
only on the deficit angle of triangles. Since the surface is embedded in R 3 , the 
numerical simulation is more time consuming than that of the model in Ref. 1281 
However, the metric is simplified by fixing the edge length L of the triangle in M 
such that L — l, and hence the numerical simulations are greatly simplified compared 
to those for the model in Ref. [55J 

We should note that the surface area of X(M) varies if the parameter a in 
S(X) = aSi+bS2 is varied while b is fixed. However, the parameter a can always 
be fixed to a — 1 because of the scale transformation X X' = (1/ y/a)X. Indeed, 
S(X) changes from S(X) = aSi+bS2 to S(X) — Si+bS2 while the partition function 
Z remains unchanged; it changes up to a multiplicative constant, under this scale 
transformation. This implies that a can always be fixed to a = l. 



2. Discrete surface model 

A discrete Hamiltonian and the partition function are introduced in this section. 
These discrete quantities are closely related to those in Ref. [28l where the two- 
dimensional surface M is embedded in R 2 by a mapping X, while in this paper M 
is embedded in R 3 . 

A two-dimensional surface M is smoothly triangulated, and a triangle A is shown 
in Fig.QJa), where I/j(i = l,2,3) denotes the edge length, and = 1, 2, 3) denotes 
the internal angle of A. The triangle Zi in Af is mapped into R 3 by X as shown in 
Fig.QJb), where the triangle X{A) is a linear one. The edge lengths and the internal 
angles of X{A) are denoted by li and (f>i, respectively. 

Let \I/ denote a coordinate mapping from A to a domain D in R 2 , then we 
have a linear triangle ^(A) in D, which is shown in Fig. QIc). The edge length 
Li(i = l, 2, 3) of A is assumed to be identical to the one in ^(A), while the internal 
angle $,(i = l,2,3) of A differs from $°(i = l,2,3) of *(Z\), where ELi $ ? = 27r - 
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(c) (d) 



Fig. 1. (a) The smooth triangle A in M, (b) the image X(A) of A by a mapping X from M to 
the external space R 3 , (c) a local coordinate domain D in R 2 and the regular triangle ^(A), and 
(d) a unit normal vector no of X(A) and those rii(i=l,2, 3) of the nearest neighbor triangles. 

In the conventional Regge calculus model, there is no difference between A in Fig. 
QJa) and W(A) in Fig. [Tfc) because <&, = 

As mentioned in the Introduction, one basic assumption is that 

L x = L 2 = L 3 = 1 (1) 

for all A. This implies that $° = 7r/3(i = 1, 2, 3). Another assumption is expressed 
as 

$! = $ 2 = $ 3 = (2) 

where $ depends on A, and $ is not always identical to 7r/3 as in the model 
of Ref. [28] We should note that "smoothly triangulated" does not always mean 
Ylj(j) ^i(i) =27r, where is an internal angle meeting at the vertex i of A in M 

Elf 

Since $ is not always given by $ = 7r/3, the deficit angle ip can be defined such 
that 

<p = 3$ - 7T, (3) 

or equivallently 

$ = $°(l + ^) = i(7r + ^). (4) 
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We should note that A is not always a linear triangle, because the edges of A are 
curved and the deficit angle (p is not always zero. On the contrary, the triangle 
X(A) is linear, because the deficit angle of X(A) is fixed to be zero. 
The surface model is defined by the discrete metric 

9ab= (p , F-cos$, \F\ < 1, (5) 

where F represents the deficit angle if. This g a b is the induced metric of the coor- 
dinate mapping '5. The inequality in Eq. (J5J) implies that ds 2 = J2ab9abdx a dXf, is 
positive definite. The area of A in M is given by 

A a = \Vl-F 2 . (6) 

The area Aa varies according to the variation of F although the edge length L is 
fixed. Since g varies as a function of F, we call the model as a Regge calculus model 
even though L is fixed. 

The Hamiltonian is defined by 

S(X,g)=S 1 + bS 2 , (7) 
Si = IJ2 S ^ A )/ A ^ S 2 = l -Y,S 2 {A)/A A , 

A A 



with 



and 



Si (A) =l(+£i + t 3 -F A (£ x £ 2 cos fa + £ 2 £ 3 cos fa + £ 3 £i cos <j> 2 ) (8) 



S 2 {A) = 3-n -ni-no-n2-no-n 3 

-F^[(ni-no)-(n 2 -no) + (n 2 -no)-(n 3 -no) + (n 3 -n )-(ni-n )], (9) 

where the symbol fa is an internal angle of X(A) as mentioned above. F A in Eqs. 
([8|) and ([9]) denotes that the value of F = cos$ depends on the triangle A. The unit 
normal vectors n^i — 0, 1, 2, 3) are shown in Fig. [Tfd). 
The partition function Z(b) is defined by 

Z(b) = Jvg J VX exp [S(X, g)] . (10) 

The integration J T>X in Z(b) is given by the multiple 3-dimensional ones such that 

VX= l[dX u (11) 
•* i=l 

where f denotes that the center of mass of the surface is fixed to the origin of R 3 . 
The symbol J T>g in Z(b) is defined by 

V 9 = I \\ dF i / II d P<- ex P (~^fS f - ^ p aS p a ~ ^pS p ) , 

i— 1 i— 1 

\Fi\<l, Pi >0. (12) 
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The integration J Y\i=i dpi is included in J Vg although p is a variable independent 
of the surface geometry. The symbol Nt = 2N—4 in Eq. (fT2|) is the total number of 
triangles. In the exponential factor in Eq. (fl"2j) . Sf is the interaction term for the 
variable F given by 

Sp^lFi-Fi], (13) 

ij 

where represents all nearest neighbor triangles i and j. In the presence of the 
interaction Sf, the metric g a b can be smoothened as a function on the surface M. 
Thus, it is physically natural to include Sf as a measure term. We should note that 
the interaction \F{ — Fj\ in Eq. (|13l) can also be expressed by (dF) 2 , however we 
use the expression \Fi — Fj\ for numerical simplicity. 
S p a and S p in Eq. (IT2t are given by 

S P A = ^2paA a , S p = ^2\pi- pj\, (14) 
A ij 

where pa (or pi) is the conjugate variable to the area Aa (or Aj) and can be 
called the "surface density". We should note that the surface density pa is not an 
external field but introduced as a variable field on the surface to see the in-plane 
surface deformation. 

The density pa is interconnected to the geometric variables X and g only 
through the interaction paAa- The term S p is introduced to define an interac- 
tion between the fileds pa ■ Since the variable p couples to A in S p a , then p becomes 
nonzero finite (or well-defined). We should note also that A a can vary without the 
density field p because Aa depends on F which varies on the surface M. However, 
the interaction S p a is expected to influence the in-plane deformation of A a ■ 

By including the measure terms of the exponential factor of Eq. (fT2f in the 
Hamiltonian, we have the effective Hamiltonian such that 

S (X, F, p) = S\ + bS 2 + XfSf + ^ p aS p a + XpSp. (15) 

Note also that the conventional Hamiltonians S\='^2 li -(Xi~-Xj) 2 and <5 f 2 = X)ij(l — 
rii • rij) are restored up to irrelevant multiplicative factors if F — 0. In this case, 
we have Sf = and A a — 1/2, and therefore we find that both S p a and S p are 
independent of the surface and can be neglected. 

3. Continuous surface model 

The surface model of Helfrich and Polyakov is defined by a mapping X from a two- 
dimensional surface M to R 3 , described by X : M 9 {x\,X2) X(x\,X2) G R 3 . 
The surface M is assumed to be of sphere topology in this paper just like in Ref. EU 
The symbol {x\, x 2 ) denotes a local coordinate of M . The image X(M) corresponds 
to a real physical membrane, however, the self-avoiding property is not assumed in 
X(M), and hence the mapping X is not always injective. 
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The Hamiltonian of the model is given by a linear combination of the Gaussian 
bond potential S\ and the extrinsic curvature energy S2 such that 

S = S 1 +bS 2 , 

51 = J Vgd 2 xg ab d a x»d b x», (16) 

52 = ~ J ^d 2 xg ab d a n»d b n», 

where b[kT] is the bending rigidity HQI11H2] The matrix <? a (,(a, b = 1, 2) in Si and 
S% is the metric on Al, g is the determinant of g a b, and g ab is its inverse. The 
symbol n M in £2 is a unit normal vector of the surface. The continuous Hamiltonian 
is invariant under the conformal transformation such as g a b — ^ g' a b = fgab for any 
positive function /, and it is also invariant under the reparametrization such as 
x — s- x' , which changes both of the variables X and g. These symmetries allows 
us to use a constant metric g a b — S a b, however, in the case of constant metric 
9ab = Sab the in-plane deformation of M is prohibited, and as a consequence the 
area of surface X(M) remains constant as long as the surface tension coefficient a in 
S = aSi+bS2 is fixed. No information of the in-plane transformation in membranes 
is obtained from such a model with g ab = S a b- 

The discrete Hamiltonians in Eqs. (0-([9]) are obtained from the continuous 
one in Eq. (fl6|) as follows. On the triangulated surfaces shown in Figs. QJa)-[lJc), 
the partial derivatives in 5i in Eq. (Q~6|) can be replaced by d\X^ — > X^ — X^, 
&2X^ — > Xtf — Xtf, where Xj 1 (g R 3 ) denotes the position of the vertex i such that 
£1 = \X%-X?\, £ 2 = \X(f-X£\. The derivatives in S 2 in Eq. (H) can also be replaced 
by din 11 — > no — 112, d 2 n^ — > no — ni, where n^(i = 0, 1,2,3) are shown in Fig.QJd). 

We make S\ and S2 to be symmetric under the permutation of the indices of £i 
and n^ such that 1 — >• 2, 2— > 3, 3 — >• 1. By including those terms which are cyclic 
under the permutation, and by multiplying a factor 1/3, we have the discrete S\ 
and S 2 in Eqs. 

This symmetrization is necessary, because if it were not for the symmetrization, 
the Hamiltonian becomes dependent on the choice of local coordinates on the trian- 
gles. Thus, the symmetrization 1 — > 2, 2 — > 3, 3 — > 1 of Si and S2 is considered to 
be a lattice analogue of the reparametrization invariance of the continuous Hamil- 
tonian in Eq. p^|) . We should note that the component F(= cos$) of g remains 
unchanged under the symmetrization. The reason of this is because the internal 
angle <3? is assumed to be independent of the vertices just as shown in Eq. @ . 

4. Monte Carlo technique 

The Metropolis Monte Carlo simulation technique is employed to update the vari- 
ables X, g and p. The update of these variables is accepted with the probability 
Min[l,exp(-*S)], where 8S=S(new)-S(old). 

The vertex position X is updated such that X —> X' = X + SX with a random 
vector SX in a small sphere. The function F is updated such that F' = F+8F, where 
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SF (e [—0.5,0.5]) is a random number, under the constraint |F'| < 1. The density 
field p is updated by p ^ p' = p + Sp, where 5p (e [—0.5, 0.5]). The constraint p' >0 
is imposed on this update. The acceptance rate for X is approximately 60%, and 
those for F and p are 80% and 90% respectively. 

One Monte Carlo sweep (MCS) consists of N sequential updates of X, Nt 
sequential updates of F, and Nt sequential updates of p. The total number of MCS 
performed after sufficiently large number of thermalization MCS is 1.5 x 10 9 ~ 
2 x 10 9 at the transition region on the TV = 8412 and N — 12252 surfaces, and 
relatively small number of MCS is performed at non-transition region and on smaller 
surfaces. 



5. Simulation results 

To see the dependence of physical quantities on the bending rigidity 6, we fix the 
parameters Xp, A p a and A p to the following values: 

(case 1), 
(case 2). 



X F =0, 
\ F =2, 



A p a = 1 
ApA = 2, A p 



A P = 1 



(17) 



In case 1, the interaction term Sf is neglected while the density term S p a and the 
interaction term S p are included in the Hamiltonian. In case 2, both terms Sf and 
S p a are included while S p is neglected. 
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Fig. 2. The mean square size X 2 vs. b in (a) case 1 and (b) case 2, the variance C X 2 vs. 6 in (c) 
case 1 and (d) case 2, and log-log plots of the peak C™J X vs. N in (e) case 1 and (f) case 2. The 
solid lines in (a), (b), (c) and (d) are drawn by the multi-histogram re- weighting technique. The 
error bars on the symbols denote the standard errors obtained by the binning analysis. 
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We firstly show in Figs.[^a) and[5fb) the mean square size 

* 2 = ^£(^-*) 2 . *=4E* ( 18 ) 

i i 

vs. b. The transition region b in case 1 is slightly smaller than that in case 2, however, 
the behavior of X 2 against the variation b in case 1 is almost identical to that in 
case 2. The error bars on X 2 is the standard errors obtained by the binning analysis, 
and the solid lines connecting the data symbols are drawn by the multi-histogram 
re- weighting technique EH. The large errors shown at the transition region in Figs. 
[2fa) and[2](b) imply that X 2 discontinuously changes. 

Figures [U[c) and[2Jd) show the variance Cx 2 of X 2 defined by 

C x , = ^{{X 2 -{X 2 )f). (19) 

We see that the peaks C™ ax increase with increasing N. To see the order of the 
transition, we show CSf* vs. N in Figs.Ete) and[2jf) in log-log scales. The straight 
lines are drawn by a least squares fitting, and we have 

C™ 2 ax ~ N a , a= 1.19 ± 0.04 (case 1), 

a =1.32 ±0.06 (case 2). (20) 

The results in Eq. (|20p indica te that a is slightly larger than 1 in both cases. From 
the finite-size scaling theory | 33 | 34 | 35 | a = \ [a < 1) implies that the transition 
is of first (second) order. Therefore, the results in Eq. (|2"0)) imply that the model 
undergoes a discontinuous transition between the smooth sphe rica l and collapsed 
phases just like the model without the metric degree of freedom EH. 

Large errors seen in 5*2 /Nb also imply a discontinuous change in S^/Nb at the 
transition region (Figs. G3a), EJb)). The specific heat is defined by 

C S2 = p(S 2 ~(S 2 )) 2 ). (21) 

We see an expected peak CE? ax in each Cs 2 and find that C^ ax increases with 
increasing N at the transition point b where Cx 2 has its peak. The C!? ax vs. N 
are plotted in Figs.[3je) and[3jf) in log-log scales. The straight lines are drawn by 
fitting the largest three data, and we have 

C™ ax ~7V ff , ct = 0.92 ±0.09 (easel), 

(7 = 1.19 ±0.08 (case 2). (22) 

We find that both of the results in Eq. (|22|) satisfy cr ~ 1 and are consistent with 

the result of the conventional model in Ref. EH This indicates that the model 

underg oes a firs t-order transition of surface fluctuations from the finite-size scaling 
tncor [331341351 

In order to see a discontinuity in S^/Nb more clearly, we show the variation 
of S2/NB against MCS in Figs. IU(a) andQJb) at the transition points. The corre- 
sponding distribution (or histogram) h(S2) of S^/iVe are shown in Figs. HJa) and 
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1.15 1.16 k 1.17 1.15 1.16 k 1.17 10 3 N 10 4 



Fig. 3. The bending energy S2/NB vs. b in (a) case 1 and (b) case 2, the specific heat C's 2 vs. b 
in (c) case 1 and (d) case 2, and log-log plots of the peak C!? ax vs. N in (e) case 1 and (f) case 2. 
The largest three data are used for the least squares fitting in both (e) and (f). 

SJb). We find a double peak structure in these ^(S^)- This structure depends on N 
and becomes more apparent as N increases, although the dependence of ^(S^) on 
TV is not shown in the figures. This implies that Sv/Nb discontinuously changes at 
the transition point and indicates that the model undergoes a first-order transition 
in both cases. We should note that the transition occurs only four times during 
1.9 x 10 9 MCS at the transition point of the N = 12252 surface (Fig. E[b)). This is 
the reason for the large errors in data X 2 and S^/iVe as mentioned above. Reliable 
estimate of critical exponents is very difficult at the first-order transition point. 

The interaction energy Sf/Nb vs. b is plotted in Figs.[5]Ja) and[5jb). The peak 
values C™ ax of the variance Cs F = (1/N){(Sf — {Sf)) ) are plotted against N in 
Figs. (He) and[UJd) in a log-log scale. The largest three data are found to scale such 
that 

C£ iax ~ N 1 , 7 = 0.18 ± 0.13 (case 1), 

7 = 0.13 ±0.08 (case 2). (23) 

These values satisfy 7 < 1 and hen ce indica te that model undergoes a continuous 
transition of in-plane deformation I 33 | 34 | 35 | ^p ne sur f ace density p, which is not 
plotted in the figures, behaves just like Sf/Nb in Figs.OJa) andE^b). However, we 
find no power low scaling behavior in C™ ax in contrast to C™ ax in Eq. ((23|) . 

The surface area defined by Eq. (|6|) is an internal geometric variable, and hence 
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it is interesting to see whether or not the phase transitions are reflected in Aa- 
Figures [6j a) and[6jb) show Aa^= YIa Aa/Nt) vs. b in case 1 and case 2. We see 
that the variation of A a against b is quite analogous to that of S^/A^b shown in Fig. 
[3l although we see that the change of A a is very small compared to the value Aa 
itself. The variance Ca a defined by Ca a — {^/N)((J2a Aa~ (J2a A a)) ), which is 
not shown as a figure, is also analogous to Cx 2 - However, the "internal" transition 
characterized by the fluctuations of Aa is not identified as a phase transition. In- 
deed, both of the data C™^ x vs. N do not satisfy the scaling relation C™^ x ~ 
Nevertheless, the surface area Aa almost discontinuously changes against the vari- 
ation of b at the transition region (Figs. EIa)|6jb)). 

We should note that the variation of the area cla(^J2a a A/Nx) of the surface 
X(M) in Figs-EJc) andEJd) is also small compared to the value of a a itself. How- 
ever, a a varies almost discontinuously just like A a - This implies that the transition 
of shape transformation is reflected in the in-plane surface deformation. 

The behavior of F shown in Fig. [6fe) is identical to that of A a in Fig. EJa), 
while the behavior of F vs. b in Fig. |6jf) is different from that of A a in Fig. [6f b) . 
We consider that the interaction described by Sf in Eq. (fl3| is an origin of this 
difference. Thus, we find that the in-plane phase structure is sensitive to the measure 
factor in Eq. (fT2]) such as the factor of Sf- 

Finally in this section, we show the Gaussian bond potential Si/N vs. b in Figs. 
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Fig. 5. The interaction energy Sp/Ng vs. b in (a) case 1 and (b) case 2, and log-log plots of 
CJJ" vs. N in (c) case 1 and (d) case 2. The largest three data in (c) and (d) are fitted to Eq. 

d2>. 




0.13 



0.128 



0.433 

A A 
0.44: 




1.04 1.05 5 1.06 1.04 



1.05 b 1.06 1.04 

F 

0.342 



1.05 b 106 



0.4415 



0.441 




case 2 



N=8412 



0.341 




1.15 1.16 J, 1.17 



1.15 1.16 J, 117 



1.15 1.16 b 1.17 



Fig. 6. The area A/± vs. b in (a) case 1 and (b) case 2, the area a/\ vs. b in (c) case 1 and (d) 
case 2, and the variance the variable F vs. b in (e) case 1 and (f) case 2. 



Ela) andlT^b). We see that the expected relation Si/N^3/2 (N— >oo) is satisfied in 
both cases. This implies that the simulations are successfully performed. We should 
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Fig. 7. The Gaussian bond potential Si/N vs. b in (a) case 1 and (b) case 2. The solid lines are 
drawn as a guide to the eyes. 



emphasize that in both cases Si is different from either A a or a a, and hence A a 
and a a are not always constrained to be a constant. This is in sharp contrast to the 
case of the conventional model, where the Gaussian bond potential J2ij(X{ — Xj) 2 
remains constant and corresponds to the surface area. 



6. Summary and Conclusion 

In this paper, using the canonical Monte Carlo simulation technique, we have studied 
a triangulated surface model, in which the metric is assumed as a dynamical variable. 
A deficit angle for triangles of triangulated surface is assumed as the metric degrees 
of freedom. We focus on whether or not the collapsing transition accompanies an 
in-plane surface deformation. It is also interesting to see whether or not the phase 
structure of the model is identical to the one of the model with a constant metric 
such as the Euclidean metric g a b — 5 a b or the induced metric g a i, = d a X fl di l X fJ ' of the 
mapping X : M — > R 3 . 

We find that the model undergoes a first-order transition between the smooth 
spherical phase at b— >oo and the collapsed phase at 6— >0. The transition is almost 
identical to the one observed in the model with constant metric in Refs. [29l |3T] 
This indicates that the transition of shape transformations is not influenced by the 
metric variable in Eq. ([5]). We also find that the phase transition accompanies a 
continuous in-plane transition, which is confirmed by the finite-size scaling analyses 
for the peak values of the variance Cs F ■ This in-plane transition is reflected in an 
almost discontinuous change in both of the internal and external quantities such 
as the surface density and the external surface area, although the discontinuous 
changes are very small compared to the values themselves. 
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